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We report the first observation of two Cabibbo-suppressed decay modes of the 5° meson. Using a 
sample of pp collisions at -^/s = 1.96 TeV corresponding to 5.9 fb~^ of integrated luminosity collected 
with the CDF II detector at the Fermilab Tevatron, we search for new B° decay modes in a sample 
of events containing J/ip tJ'^l^~ decays. We reconstruct a B° — > J/i/j K*{8Q2)'' signal with 
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K*{892)^ — >■ i('+7r~, observing a yield of 151 ± 25 events with a statistical significance of 8. Oct. 
We also reconstruct a 5° J/ip K% signal with K% — > 7r"'"7r~, observing a yield of 64 ± 14 events 
with a statistical significance of 7.2ct. From these yields, we extract the branching ratios B{B1 — >■ 
J/ipK*{i92f) = (8.3 ± 3.8) X 10"^ and B{B° Jj'i'K^) = (3.5 ± 0.8) x 10"^ where statistical, 
systematic, and fragmentation-fraction uncertainties are included in the combined uncertainty. 

PACS numbers; 14.40.Nd, 12.15.Ff, 12.15.Hh, 13.20.He 



I. INTRODUCTION 

This paper presents the first observation of the 
Cabibbo-suppressed decays -> J/V'-ft:*°(892) 
and — )■ J/V' K'g (and the corresponding 
charge conjugate decays) using a sample de- 
rived from an integrated luminosity of 5.9 fb^^ 
of proton-antiproton collisions at sfs — 1.96 TcV 
produced at the Fcrmilab Tevatron. In addi- 
tion to isolating these signals, we normalize the 
observed yields to the corresponding Cabibbo- 
favored i?° decay modes (S" J/ipK*'^ , where 
K*^ refers to iv:*°(892), and B° J/^K'^) to 
extract the branching ratios for these newly ob- 
served Bg decay modes using the relation 



g(i?g ^ J/^K) 
B{B" ^ J/^K) 



(1) 



where K represents Kg or X*", Arei is the rel- 
ative acceptance, fs/ fd is the ratio of fragmen- 
tation fractions and N{B° J/^pK)/N{B" -> 
J/il;K) is the measured ratio of yields. 

In the naive spectator model, the ratio 
of branching ratios is given by the ratio of 
the squares of the Cabibbo-Kobayashi-Maskawa 
(CKM) elements 



cd\ 



B{BO^J/iPK) iKsp 



= 0.051 ± 0.006, (2) 



which is derived from | V^d] — 0.230 ± 0.011 and 
\Vcs\ = 1.023 ± 0.036 Q. Assuming a rela- 
tive acceptance Arei of unity, we estimate a ra- 
tio of yields. The value for fs/fd is extracted 
from the most recent CDF measurement [1] of 
fs/ifu + fd) X B{Ds -)> (/)7r) and fu/fd along with 
the current world-average value [i| for B{Ds — 
(j)Tr). Combining the value fs/fd = 0.269 ± 0.033 
with Eq. [2] yields 
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iV(gg ^ J/i:K) ^ BjB'i^ J/^K) fs 1 
N{B° ^ J/ijK) B(B° ^ J/ijK) fd Arei 
= 0.014± 0.002. 

(3) 

While the result holds only in the simple spec- 
tator case, it provides useful guidance that we 
might expect one to two Cabibbo-suppressed 
S° J/i)K events for every 100 Cabibbo- 
favored B^ J/ipK events. 

With the establishment of the decay modes 
presented here, future measurements can be con- 
sidered that will further aid our experimental in- 
vestigation into the physics of the 5° system. 
The success of the CKM three-generation de- 
scription of charge conjugation-parity inversion 
( CP) violation Q in the bottom and kaon sectors 
has continued to motivate additional, more pre- 
cise tests of CP violation in the flavor sector. In 
recent years, attention has turned to the B^ me- 
son as new territory to explore the possibility of 
non-standard-model contributions, specifically in 
the CKM matrix element Vts- Precise measure- 
ment of the frequency of B^ flavor oscillations [3| 
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has significantly limited the magnitude of new 
physics amplitudes. However, possible large new 
physics phases remain poorly constrained. 

Cabibbo-suppressed i?^ modes could provide 
complementary information on the mixing 
phase and on the width difference AF^o = 
Fgo^ — r^o^ where r^o^(r^o^) is the width 
of the light, even (heavy, odd) CP eigen- 
state Q. The decay 5° ^ J/V^K*(892)° is a 
pseudoscalar to vector- vector transition and can 
be used to help disentangle penguin contributions 
in — 7- J/^p (j) . With a sufficiently large data 
sample, it would be possible to measure Ar^o 
and the polarization amplitudes. Furthermore, 
the Cabibbo-suppressed decay B^ — > J/ij^Kg 
is a CP-odd final state (ignoring CP violation 
in the kaon system) and therefore a measure- 
ment of the lifetime in this decay mode is a di- 
rect measure of F^o = 1 /r^o . With a larger 
data sample, a tagged CP asymmetry analy- 
sis of the B^ — >■ J/ipK'g mode, in conjunction 
with our precise knowledge of CP violation in 
_B" J/ipKg, can yield information on the an- 
gle 7 of the unitarity triangle ^ . 

After a description of the detector, data sam- 
ple, and simulated samples utilized here, we de- 
scribe the B° J/ipK*°{892) analysis in Sec. 

III, followed by the B'^ — J/ipKg analysis in Sec. 

IV. Section V then describes the acceptance cal- 
culation for both modes, followed by the results 
in Sec. VI. 



II. CDF DETECTOR, DATA, AND 
MONTE CARLO SAMPLES 

The data used in these analyses correspond to 
an integrated luminosity of 5.9 fb~^ and were col- 
lected by the CDF II detector from March 2002 to 
February 2010 using di-muon triggers. The CDF 
II detector is a general purpose, cylindrically 
symmetric detector. A more detailed description 
can be found elsewhere Q- The sub-detectors 
relevant for these analyses are briefly discussed 
here. Charged particle trajectories (tracks) are 
measured by a system comprised of eight layers 
of silicon microstrip detector (SVX) and an open- 
cell wire drift chamber (COT), both immersed in 
a 1.4T axial magnetic field. The silicon detec- 
tor Q extends from a radius of 1.5 cm to 22 cm 
and has a single-hit resolution of approximately 
15 fim. The COT drift chamber 9] provides up 
to 96 measurements from radii of 40 cm to 137 cm 
and covers the range |?7| <1 [l3]. Combined 



COT-I-SVX charged particle momentum resolu- 
tion is cTp^/pT^ = 0.07%[GeV/c]-i. Outside 
the calorimeters reside four layers of planar drift 
chambers [ll| (CMU) that detect muons with 
transverse momentum pT >1.4GeV/c within 
\r]\ < 0.6. Additional chambers and scintilla- 
tors [13 (CMX) cover 0.6< \ri\ <1.0 for muons 
withpT >2GeV/c. 

The di-muon triggers collect a sample of 
J/i(j fJ-^IJ'~ candidates. At the first level 
of a three-level trigger system, an electronic 
track processor (XFT) [13| uses COT informa- 
tion to find tracks and extrapolate lj| those 
with pt >1.5(2.0) GeV/c to track segments in 
the CMU (CMX) muon-chambers. Events pass 
this first trigger level if two or more XFT tracks 
arc matched to muon-chamber track segments. 
The second trigger level requires those tracks to 
have opposite charge and an appropriate open- 
ing angle in the plane transverse to the beam- 
line. Finally, at level 3, full tracking information 
is used to reconstruct J/ip ^ A*"*"/^" candidates. 
Events with a candidate in the mass range 2.7 to 
4.0GeV/c^ are accepted. 

To identify B° and 5^ decay candidates, we 
pair J/ip candidates with Kg — > tt+tt" and 
j^*o _^ iir+7r~candidates. The reconstruction 
of Kg — > TT+TT^ and if*" K+tt^ candidates 
starts from pairs of oppositely-charged tracks 
fit to a common interaction point (vertex). In 
the — >■ J/ip Kg analysis, we reconstruct two 
tracks as pions and combine them to define a Kg 
candidate, where the invariant mass of the two pi- 
ons is constrained to the known Kg mass Ij . In 
the — > J/ipK*'^ analysis, we reconstruct the 
if *° candidate from the combination of a tt and a 
K. If two A' *° candidates are reconstructed with 
the same tracks, with the only difference that the 
kaon and pion hypotheses are interchanged, we 
select the K*° candidate whose mass is closest 
to the pole value of 896MeV/c^. We perform a 
kinematic fit of each B candidate where the fi- 
nal state tracks are constrained to come from a 
common decay point and the invariant mass of 
the muon pair is constrained to the known J/ip 
mass |l| . These preliminary selection criteria for 
B° and i?^ candidates are listed in Table H] Ad- 
ditional selection criteria optimized for the indi- 
vidual channels are described in Sees. Ill and 
IV. 

Simulated samples of i?° and B^ decays are 
used to optimize event selection, model signal 
distributions, and assess systematic uncertain- 
ties. For our default Monte Carlo simulation 
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(MC) samples, we generate single b hadrons ac- 
cording to the predicted next-to-leading order 
QCD calculation [l^. For systematic studies, 
we also generate single b hadrons according to 
momentum and rapidity spectra measured by 
CDF These hadrons are then decayed us- 
ing the EVTGEN package fl^ and then fed into 
a GEANT simulation of the CDF detector [l7| . 
The simulated data are then processed and re- 
constructed in the same manner as the detector 
data. In the case of J/^K*° mode, it IS nec- 
essary to specify the polarization parameters in 
the simulation. For both i?° and we use 
transversity basis [111 polarization amplitudes 
1^0 P = 0.6 and \A±\'^ = 0.22, which are simi- 
lar to the PDG values of \An^ = 0.571 ± 0.008 
and \A^\'^ = 0.22 ± 0.013 Q. For systematic 
acceptance studies, MC samples with other po- 
larization values were generated. 

In all of the MC samples generated, and 
throughout the analyses presented below, we as- 
sume that there is no CP violation in B'^ mix- 
ing or decay. We additionally assume that equal 
numbers of and mesons, as well as equal 
numbers of 5° and B^ mesons, are produced in 
the pp collisions. From these assumptions, this 
untagged analysis is insensitive to CP violation 

decays, and the width difference in the 
system is given by AF go . 

III. B° ^ J/ipK*° ANALYSIS 

We optimize the selection criteria to provide 
the highest likelihood for evidence of this mode. 
This is done by maximizing S/ (1.5 -I- Vb), where 
S refers to the number of signal events and B 
is the number of background events in the sig- 
nal region. Reference jl^ demonstrates that 
this quantity is well suited for discovery. For 
the signal sample, a B° — >• J/ipK*'^ MC sam- 
ple is used. For the background sample, we 
use J/tpK*'^ candidate events from data with 
the requirement that the reconstructed candidate 
mass Mb falls in the range 5.6GeV/c^ < AIb < 
5.8GeV/c^. This "upper sideband" region con- 
tains events kinematically similar to the com- 
binatorial background in the signal region and 
is not contaminated by residual signal events. 
We optimize simultaneously over the transverse 
momenta pt{t^^) and pt{K^), the B^ trans- 
verse decay length Lxy{B'^), and the 5° decay 
kinematic-fit probability. The final cuts we use 
arepT(7r~) > 1.5 GeV/c, pt{K+) > 1.5 GeV/c, 



Lxy{B^) > 300 /im and fit probability greater 
than 10"^. 

Particle identification using specific ionization 
[dE/dx) in the COT was evaluated to further 
separate K*^ — >■ K^tt~ from tt+tt" and K^K~ 
backgrounds. Although further background re- 
duction could be achieved, the corresponding 
reduction in signal efficiency rendered particle 
identification unprofitable, and we choose not to 
use it. 

We determine the B'^ and i?° yields using a 
binned likelihood fit in the candidate masses. We 
model the signal contributions with templates 
composed of three Gaussians obtained from fits 
to B'^ MC. The two dominant, narrow Gaus- 
sians model detector resolution effects and also 
account for cases where the identities of the tt 
and K from the decay are interchanged. As 
mentioned above, we identify events where both 
TT-K and K-tt hypotheses pass the selection cri- 
teria and, in those cases, choose the combination 
closest to the A'* (892) mass to ensure that candi- 
dates are not used twice. Approximately 10% of 
B -J> J/ipK*° events are reconstructed with the 
incorrect tt-K assignment. These events peak at 
the B masses, but have a significantly broader 
width. 

A wide Gaussian models misreconstructed sig- 
nal events and other non-Gaussian resolution ef- 
fects. The relative contributions, means, and 
widths of each Gaussian are fixed in the fit. The 
5° templates used in the fit are identical to B^ 
templates, except for a shift of 86.8 MeV/c^ in 
the mean value of the three Gaussians. This 
value corresponds to the known [T, mass dif- 
ference between i?° and The MC slightly un- 
derestimates the mass resolution, so the widths 
of the two narrow Gaussians are multiplied by 
a scale factor common to the i?° and B'^ tem- 
plates, which is allowed to float in the fit. The 
scale factor is not applied to the third Gaussian 
since it is not expected to be affected by detector 
resolution effects as the other two are. Moreover, 
a common mass shift is added to the means of 
all Gaussian templates to account for a possible 
mass mismodeling in the MC. This mass shift is 
floating in the fit. 

The J/ijj K*^ analysis has three primary 

background contributions: events with random 
track combinations (combinatorics) , partially re- 
constructed b hadrons, and B^ — > J/ipcj) decays. 
Combinatorial background arises from sources 
such as a real J plus two other tracks, where 
the J/ij} could be either prompt or coming from 
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a B decay. Another source arises from false 
J/ ip candidates reconstructed from misidentified 
hadrons. The combinatorial background is mod- 
eled in the fit with an exponential function. 

Backgrounds from partially reconstructed b 
hadrons come from multibody decays where a tt, 
K, or 7 is not reconstructed, for example, the 
decay mode B° J/ipK*°n°. We fit this back- 



ground with two ARGUS functions [2l|, one for 
partially reconstructed and another for par- 
tially reconstructed B^. The ARGUS function 
parameterization for m < mo is 



/(m) = TVi X 



X e 



(4) 



where mo is the mass cutoff, C the decay con- 
stant, and A'^i is the normalization. The function 
is zero for m > mo. The ARGUS function for 
partially reconstructed 5° has a fixed mass cutoff 
of m{B°) - m{TT°) = 5.140 GeV/c2 and the func- 
tion for partially reconstructed B^ has a fixed 
mass cutoff of m(B°) - m{n°) = 5.220 GeV/c^. 
The decay constants of the two functions are 
constrained to be the same, and the normaliza- 
tions are independent. Each ARGUS function 
is convoluted with a Gaussian having a width of 
12MeV/c^ to account for detector resolution ef- 
fects. 

Since it is possible for — > J/i' 4> candi- 
dates to pass the J/ipK*^ reconstruction crite- 
ria, B^ — J/ip(j) must be considered as a back- 
ground. We use a template consisting of two 
Gaussians, extracted from simulation, to model 
this background in the J/ipK*'^ fit, where both 
Gaussians are primarily modeling detector res- 
olution effects. We fix the widths, means, and 
relative contributions from each Gaussian in the 
final fit. We multiply the constant width of 
the narrower Gaussian by the same scale factor 
used in the signal templates. We constrain the 
J/ij) (j) contribution in the J/t/j K*° fit by 
measuring the yield of i?° — >■ J/ipcj) in the data 
using selection criteria efficient for reconstructing 
B^ J/ip(f>. We then use simulation to calcu- 
late the fraction of those J/tp(j> events that would 
satisfy the J/'ipK*^ selection. 

We perform a binned log likelihood fit to the 
J/ip Ktt invariant mass distribution using the 
templates for signals and the background func- 
tions described above. The mass distributions in 
data for J/ip K*^ candidates and the final fit ap- 
pear in Fig. [H The yields for B" J/%IjK*° 



and B° J/i}}K*° signal are 9530 ± 110 and 
151 ± 25 respectively. The ratio iV(S^ -J> 
J/ijK*°)/N{B° ~> J/i;K*°) is 0.0159 ± 0.0022 
(stat). 

We determine the statistical significance of the 
-> J/^K*° signal by fitting the mass distri- 
bution without the contribution (background- 
only hypothesis). For likelihood C, we inter- 
pret — 21og£ as a distribution. We use Ax^ 
with one degree of freedom to determine that the 
probability of background fluctuations produc- 
ing a comparable or greater signal is 8.9xl0~^^ 
or 8.0(7. This is the first observation of the 
J/i}jK*° decay 

We consider several sources of systematic un- 
certainty in the measured ratio of N{B'^ — >■ 
J/iIjK*°)/N{B° ^ J/ijK*"). The modeling of 
the B^ and B^ signal peaks can infiuence the 
measurement of the ratio. To quantify the effect 
of the mismodeling we repeat the fit using two 
Gaussian templates instead of three for the sig- 
nal. The fit value of N{B^^ )/N{B") is shifted by 
7x10-^^. 

We vary the input mass difference between 
B'^ and i?s in the templates within its uncer- 
tainty of 0.7 MeV/c^. The difference in N{B° ^■ 
J/^K*")/N{B" J/ijK*^) with the alternate 
templates is 2x10^'"'. This is sufficiently small 
that we ascribe no systematic uncertainty for this 
effect. 

The shape of the combinatorial background is 
another source of systematic uncertainty. In this 
case, we use a power function instead of an ex- 
ponential. We assign an additional systematic 
uncertainty of 2x10""' to account for this effect. 

In the likelihood fit, we allow the combinatorial 
background contribution to fioat. We performed 
a study to evaluate how the ratio of yields de- 
pends upon the specific, arbitrary choice of the 
fit range. We compare the main fit, which al- 
lows the combinatorial background to float over 
the entire fit range, to a control case where the 
combinatorial contribution is fitted in the up- 
per sideband and extrapolated to the full mass 
range prior to the final fit. Due to the differ- 
ence in the result from these two methods, we 
include a systematic uncertainty of 0.0050 on the 
iV(B^ ^ J/i:K*°)/N{B" ^ J/iPK*") ratio. 

To study the uncertainty in the B'^ J/ip (p 
contribution, we repeat the fit while doubling the 
fraction of -B^ — ?> J/4'4> candidates. The result- 
ing shift of 2 X 10"'' is assigned as the uncertainty 
in the i?° J /^jj (p contribution. 

We add the different systematic uncertainty 
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contributions, summarized in TableHIl in quadra- 
ture resulting in a final value oi N{B'^ — > 
J/^PK*°)/N{B° ^ J/'iPK*°) of 0.0159 ± 
0.0022 (stat) ± 0.0050 (syst). 

IV. B° ^ J/iP K*-^ ANALYSIS 

The i?" — > J/tp Kg decay has several differ- 
ences compared to the — J/ipK*^ decay. It 
contains a Kg, which has a relatively long life- 
time of cr=2.68 cm. We use the displacement be- 
tween the reconstructed Kg decay point and the 
reconstructed B decay point in the event selec- 
tion to reduce backgrounds such as B^ ^ J/ip cj). 
Finally, as in the B^ system, we expect the 
B^ — J/ipKg signal to be smaller than that 
of the B° J/ipK*^ mode. Therefore, we 
use a Neural Network (NN) technique to take 
full advantage of all the kinematic variables and 
their correlations. We use the NeuroBayes [l^ 
NN package. The NN provides an output value 
close to +1 for signal-like events and near —1 for 
background-like events. 

We train the NN using simulated B° MC 
events as a signal sample. We use data from the 
upper sideband in the J / Kg candidate mass 
distribution, well separated from the signal re- 
gion, as a background training sample. We use 
as inputs for the NN the quantities listed in Ta- 
ble mil These input quantities are chosen as 
variables with good discriminating power which, 
alone or in combination, do not bias the mass 
spectrum. After the training, the NN achieves 
strong discrimination between signal and back- 
ground as shown in Fig. [2^. 

As in the B° J/ipK*° analysis, we opti- 
mize the selection by maximizing S'/(1.5 -I- ^/B). 
The signal S is modeled using MC events in 
the reconstructed mass range 5.350 GeV/c^ < 
Mb < 5.400 GeV/c^. The background B is mod- 
eled using J/ Kg candidates in data populating 
the mass range 5.430 GeV/c^ < Mb < 5.480 
GeV/c^. The figure of merit suggests a cut value 
in the NN response of 0.88 as shown in Fig. ^p. 

The fitting technique is similar to the 
J/ip K*° analysis. We obtain the yields of B" — > 
J/ip Kg and — J/ip Kg signals in a binned 
likelihood fit to the invariant mass distribution. 
We again model the B'^ and 5° signal contri- 
butions with three Gaussian templates obtained 
from fitting B^ ^ J /ip Kg MC and use the mass 
difference between 5° and B'~' for the formation 
of the J/ijj K*^ template. The two major 



sources of background in this analysis are combi- 
natorial background and partially reconstructed 
fe-hadron decays. We model these with the same 
functional forms used in the i?° — J/ip K*^ anal- 
ysis. However, we include only one ARGUS func- 
tion because the contribution of partially recon- 
structed i?s is negligible. An additional back- 
ground in this analysis is Ab" — > J/V'A decays 
where the p from the A decay is assumed to be a 
TT. In order to suppress the Aj,*^ contribution, we 
apply a cut to the angular variable cos(0^o .n.^), 
where Oj^o is the angle between the Kg candi- 
date pt in the lab frame and the lower px pion 
(772) in the Kg center-of-mass frame. Cutting out 
events with cos(6';fO jr2)< -0.75 removes 99.8% of 

the Ah° while retaining 86% of the B°. The resid- 
ual Ab* contamination is less than one event and 
is neglected. The invariant mass distribution for 
J/tp Kg and the fit result including the different 
contributions are shown in Fig. |3l 

We determine the yields of the B'^ ^ J/tp Kg 
and B° J/i^Kl signal to be 5954 ± 79 
and 64 ± 14 respectively. As with the i?" — > 
J/ipK*^ case, we determine the statistical sig- 
nificance of the i?j — J/i) K^ signal by fitting 
the mass distribution without the contribu- 
tion (background-only hypothesis), a difference 
of one degree of freedom between the two hy- 
potheses. For likelihood C we interpret — 21og£ 
as a &iid use the difference in that quan- 
tity to determine that the probability of back- 
ground fiuctuations producing a comparable or 
greater signal is 3.9x10""'^'^ or 7.2a. The value of 
N{B° J/^K'^g)/N{B^ Jl-^Kl) is 0.0108 
± 0.0019 (stat). 

The sources of systematic uncertainty are sim- 
ilar to the other analysis. In this case the ab- 
solute uncertainties for the ratio are 6x10"'' 
from the combinatorial background contribu- 
tion, 6x10"'* from the combinatorial background 
modeling, 5x10""* from the signal modeling and 
1.3x10"''^ from the mass difference between i?° 
and -B". The systematic uncertainties are sum- 
marized in Table |IT1 We sum the contribu- 
tions in quadrature resulting in a total system- 
atic uncertainty of ±0.0010. The final value of 
N{B° J/iPK%)/N{B° J/^K'^) is 0.0108 
± 0.0019 (stat) ± 0.0010 (syst). 

V. ACCEPTANCE CALCULATION 

To determine the ratio of branching ratios 
B{B"^ J/tpK)/B{B° J/i^K), where K rep- 
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resents Kg or K*^, the relative acceptances of 
B° J/ipK°s to B° ^ J/V-^s and 5° ^ 
J/^K*° to B° J/^K*° need to be deter- 
mined. We use MC samples to extract Arei as 
follows: 



■^rel — 



N{B" J/il)K pass) /N{B" J/t^K gen) 
N{B^, J/i;K pass)/ NiB'^ -> J/^PK gen) 

(5) 

where iV(gen) is the number of MC generated 
signal events, -/V(pass) is the number of events 
passing all selection requirements, and K repre- 
sents K°s or if*". 

We determine the value for Arei to be 1.057 
± 0.010 for the K*° channel and 1.012 ± 0.010 
for the Kg channel. We determine the statisti- 
cal uncertainty on the acceptances for B'^ and 
B^, assuming binomial statistics. This MC sta- 
tistical uncertainty is reported as a systematic 
uncertainty on Arei- 

The data sample utilized in this analysis was 
acquired using a number of variations on the 
J ftp — T' ^J■'^^^~ trigger. We have verified that the 
acceptance calculation is robust and consistent 
across all kinematic variations of these triggers. 

Several other effects contribute to the system- 
atic uncertainty on Arei- Uncertainty in i?^ and 
B^ lifetimes introduce an uncertainty on the ac- 
ceptance through the transverse decay length re- 
quirement. For i?° — >■ J/ijjKg analysis, we gen- 
erate different MC samples, varying the lifetimes 
by one standard deviation with respect to their 
measured values. We use the average measured 
value for B^ and the evaluated r^o value for 
Bg Q. The maximum deviation of Arei is 0.028, 
and we take this value as a systematic uncer- 
tainty. 

For the B° J/'iJjK*° analysis, the proce- 
dure to evaluate the systematic uncertainty is 
slightly different. The B° J/ipK*" decay is 
an unknown admixture of CP-even and CF-odd 
states which have different lifetimes. The world- 
average currently gives AFbo/Fbo = 0.092l:o °54 



maximum deviation on Arei of 0.046. 

Another source of systematic uncertainty 
arises from the momentum spectra of the B° and 
B^- Since we normalize our i?° signal to the B^ 
mode, we are sensitive only to mismodeling in the 
ratio of pt(5°) versus pt{B^)^ which should be 
quite small. We compare the default i?° and i?° 
samples which use a next-to-leading order QCD 
'calculation [l^ to the pr spectrum measured by 
CDF 0. In the B° J/iPK*° analysis, the 
value of Arei varies by 0.029 when using these 
alternative production spectra and we take this 
value as a systematic uncertainty. Likewise, for 
the i?g — )• J/ip Kg analysis, the change in Arei is 
0.032. 

Our relative acceptance is calculated assuming 
that the polarization in B^ — >■ J/ipK* is identi- 
cal to the polarization in B^ — >■ J/ipK*. Since we 
have no a priori knowledge of the actual polar- 
ization in the i?" mode, we compute the system- 
atic uncertainty by allowing all possible values 
for the polarization. We generated MC samples 
1, A\\ = 1, and A± = 1. The maximum 



for A 







variation from any of these polarizations leads to 
a systematic uncertainty on Arei of 0.261. Since 
the angular distributions arising from polariza- 
tion are clearly the dominant systematic uncer- 
tainty, we have neglected the correlation between 
polarization and lifetime in assessing the uncer- 
tainties. 

Table IIVI shows a summary of the system- 
atic uncertainties on Arei for both measurements. 
Summing these contributions in quadrature, we 
find Arei = F057 ± 0.010 (stat) ± 0.267 (syst) 
for the analysis and Arei = 1-012 ± 0.010 
(stat) ± 0.042 (syst) for the K^ analysis. 

VI. RESULTS 

Using the values of Arei described above, we 
find 



for F^o = ki^B" + Tb" ) where F^o and 



/,S(i?0-> J/Vif*°) 



3? = BO -r 1 BOJ UUJ, wneie 1 go^ anu ^^^(^o ^ J/^R*") 

are the widths of the heavy and light mass 



0.0168 ± 0.0024(stat) 



eigenstates respectively. If the were either 
all B'^j^ or B'^j^, the maximum lifetime change 
would be 5%. To evaluate the effect on Arei, 
we reweig ht the default B° ^ J/^PK*° lifetime 
distribution. The reweighting is performed by 
normalizing the default lifetime distribution and 
comparing it to distributions with the lifetime 
increased or decreased by 5%. This leads to a 



± 0.0068(syst) 



(6) 



and 

/.g(i3°^ J/^gg) 



/.mi?o^^/^;,0) - 0.0109 ±0.0019(stat) 



±0.0011(syst). 



(7) 
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To determine the ratio of branching ratios, 
we combine these results with the most recent 
CDF measurement of fs/{fu + fd) x B{Ds 
4>tt) and fu/ fd with the current world- average 
value [l for b{d, <l)Tr) to yield f^/fd = 0.269 
± 0.033. We quote the systematic uncertainty 
coming from the fs/fd uncertainty as "frag". 
The ratio of branching fractions to the reference 
decays are: 



In conclusion, we present the first observa- 
tion and branching ratio measurement of the 
Cabibbo suppressed decays — > J/ipK*'^ and 
B^ —J' J/ipKg. With larger data samples and 
additional analysis, these modes can be used to 
further explore the properties of the i?^ system. 



0.062 ±0.009(stat) 

± 0.025(syst) ± 0.008(frag) 

(8) 



and 

B{B°^J/^K^s) 
B{B" J/ipKl) 



0.041 ± 0.007(stat) 

± 0.004(syst) ± 0.005(frag). 

(9) 



The relative branching ratios observed for both 
modes are in good agreement with the expecta- 
tion based upon the pure spectator model. 

We use the world-average values for B{B^ 
J/iIjK*°) and B{B° J/i^K^) [H for normal- 
ization to calculate the absolute branching frac- 
tions: 



B{B° J/4> K*") = (8.3 ± 1.2(stat) ± 3.4(syst) 

± l.O(frag) ± 0.4(norm)) x 10"^ 

(10) 

and 

B{B"^ ^ J/tp K°) = (3.5 ± 0.6(stat) ± 0.4(syst) 

± 0.4(frag) ± O.l(norm)) x 10"^ 

(11) 
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TABLE I: Selection criteria for B" J/i/jK candidates and 3° J/tpK candidates, where K represents 
K*° or K°s. 



Variable (Units) 


B^-^JI^IjK*'' 


B'^^^J/^K'^s 


B"/ Bs candidate four-track fit x'' 


< 50 


- 


B" / B" candidate four-track fit probability 


- 


> 10"^ 


B" / B" candidate transverse momentum pr (GcV/c) 


> 6 


> 4 


B" / B" candidate impact parameter (/im) 


< 50 


- 


5°/ B^ candidate transverse decay length significance L^y/cr 




> 2 


J/ij) candidate mass (GeV/c^) 


> 3.05 


> 2.8 




< 3.f5 


< 3.3 


J/tj) candidate 3D two-track fit 


< 30 


< 30 


K candidate mass (GeV/c^) 


> 0.55 


> 0.55 




< 0.846 


< 0.846 


K candidate 3D two-track fit 


< 30 


< 20 


K candidate transverse decay length L^y (cm) 




> 0.5 


transverse momentum pT (GeV/c) 


> f.5 


> 1.5 


A(/> between the two muons (radians) 


< 2.25 


< 2.25 


/xi charge x jj,2 charge 


= -1 


= -1 


A2; in the beam line between the two /it (cm) 


< 5 


< 5 


7r transverse momentum pT (GeV/c) 




> 0.5 



TABLE II: Systematic uncertainties for the ratio of yields. All rmmbers in percent. 



Source 


j]V(s^^j/^K*") ^o/^ 

O N(B<'^J/4,K'<i) (./"i 


rN(B",^J/'^K'-^) ™n 
" N(B0^J/^K%) \'^> 


Signal Modeling 


4.4 


4.6 


Mass difference between B^ and _B° 


0.1 


0.1 


Combinatorial background modeling 


1.3 


5.6 


Combinatorial background contribution 


31.4 


5.6 


B'^ — > J/ih d) contribution 


1.3 




Total 


3i.8 


9.2 



TABLE III: Variables used as input in the NN training. 

Input variables in the NN 
B'^/Bs candidate transverse momentum 
B'^ / Bs candidate four-track decay point fit 
B'-' /B" candidate proper decay length 
B^ /B1 candidate impact parameter 
J/V' candidate transverse momentum 
J/tp candidate mass 
J/^/i candidate proper decay length 
J candidate impact parameter 
K% candidate transverse momentum 
K% candidate mass 
K% candidate proper decay length 
K% candidate impact parameter 
TT transverse momentum 
TT impact parameter 
H transverse momentum 
/X impact parameter 

H cosine of the helicity angle in J/tj; rest frame 
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FIG. 1: (a) Invariant mass distribution in data for J/tp K*^ candidates and fit including the different contri- 
butions, (b) We enlarge the distribution in the signal region for more detail. 



TABLE IV: Systematic uncertainties for the relative acceptances. All numbers listed in percent. 



Source SArei{B^ 


J/V'A'*" (%)) 


J/7/;ifg (%)) 


Lifetime for and B'^ 


4.4 


2.8 


B hadron pr spectrum 


2.7 


3.2 


Polarization 


24.7 




Total 


25.3 


4.2 
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FIG. 2: (a) NN response where the sohd line is signal simulation and the dashed one is sideband data, (b) 
Figure of merit 5/(1.5 + \/!B) as a function of NN response. The vertical line indicates the optimized cut in 
the NN response. 
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FIG. 3: (a) Invariant mass distribution in data for J/ip K% candidates and fit including tfie different contri- 
butions, (b) We enlarge the distribution in the signal region for more detail. 



